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Abstract. We investigate the production of hadrons in nuclear collisions within the 
framework of the thermal (or statistical hadronization) model. We discuss both the 
ligh-quark hadrons as well as charmonium and provide predictions for the LHC energy. 
Even as its exact magnitude is dependent on the charm production cross section, not 
yet measured in Pb-Pb collisions, we can confidently predict that at the LHC the 
nuclear modification factor of charmonium as a function of centrality is larger than 
that observed at RHIC and compare the experimental results to these predictions. 

One of the major goals of ultrarelativistic nuclear collision studies is to obtain 
information on the QCD phase diagram pQ. The experimental approach is the 
investigation of hadron production at chemical freeze-out [2] . Hadron yields measured 
in central heavy ion collisions from AGS up to RHIC energies can be described very 
well (see [3] and refs. therein) within a hadro-chemical equilibrium model. In our 
approach the only parameters are the chemical freeze-out temperature T, the baryo- 
chemical potential //& and the fireball volume V. The main result of these investigations 
was that the extracted temperature values rise rather sharply from low energies on 
towards y / ^v r /v —10 GeV and reach afterwards constant values near T=160 MeV, while 
the baryochemical potential decreases smoothly as a function of energy. The limiting 
temperature (!] behavior suggests a connection to the phase boundary and it was, indeed, 
argued [5] that the quark-hadron phase transition drives the equilibration dynamically, 
at least for SPS energies and above. The conjecture of the tricritical point [6] was put 
forward for the lower energies. Alternative conjectures are that the thermodynamical 
state is a generic fingerprint of hadronization IB] or is a feature of the excited QCD 
vacuum [9]. 

The values of T and \x\, obtained from fits can be parametrized as a function of ^/sNN 
with the expressions: T = T Km /[l + exp(2.60 — ln(^/sAT7v(GeV))/0.45)] and /^[MeV] = 

1303/[1 + 0.286 v /s JViV (GeV)], with the "limiting" temperature T Hm =lU MeV. The fi b 
value expected at the LHC is around 1 MeV. We employ these parametrizations to 
compare the model to data over a broad energy range. As an illustration, the production 
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Figure 1. Energy dependence of the relative production ratios of protons (p) and kaons 
(K~) to pions (tv~). The open symbols represent the data points for the NA44 and 
PHENIX experiments at SPS and RHIC, respectively (the full symbols are for NA49 
and STAR data). Note that ad-hoc feed-down subtraction was applied for STAR data 
for protons (25%) and for PHENIX data for pions (10%). The errors of the p/ir and 
p/n data are smaller than the symbols. 



yields of protons and kaons relative to pions are shown in Fig. [TJ demonstrating that 
the model describes the data well (although smaller p/n + and p/ir~~ ratios are measured 
by PHENIX and at SPS [3]). The trends seen in the p/ir + and p/ir~ ratios reflect 
both the strong increase followed by saturation for T and the strong decrease of fib as a 
function of a/saT/v. While the K~ /ir~ ratio shows a monotonic increase and saturation 
as a function of energy, the K + /tt + ratio shows a maximum at a beam energy of 30 
AGeV. In the thermal model this maximum occurs naturally at as an effect of the steep 
rise and saturation of T and the strong monotonous decrease in ^ [3]. 
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Figure 2. Model prediction for hadron yields relative to pions at LHC (left panel) 
and of anti-matter to matter production at RHIC and LHC. 
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In Fig. [2] we provide predictions for the production of various hadrons relative to 
pions, to be tested soon by experiment. Preliminary ALICE data [10] indicate a lower 
p/n~ ratio. We show also how dramatic the balance between matter and anti- matter 
production is changing from RHIC to LHC energies, as illustrated by the ratios of 
anti-baryons to baryons. Prediction for (anti-)hyper-nuclei are also available [TT] . 

We now turn to the heavy-quark sector and compare the model predictions to 
data on charmonium prodution. Charmonium is considered, since the original proposal 
about its suppression in a Quark-Gluon Plasma (QGP) p2], an important probe of the 
energy density reached in the deconfined fireball produced in ultra-relativistic nucleus- 
nucleus collisions (see [13]). Because the large mass of charm quarks, heavy flavor 
hadron production cannot be described in a purely thermal approach as that discussed 
above. It was realized in [14J that charmonium production can be well described in 
the statistical model by assuming that all charm quarks are produced in initial, hard 
collisions while charmed hadron and charmonium production takes place exclusively at 
the phase boundary with statistical weights calculated in a thermal approach (for a 
recent review see [H]). An important element is thermal equilibration, at least near the 
critical temperature, T c , which we believe can be achieved efficiently for charm only in 
the QGP. In recent publications [16] we have demonstrated that the data on J/^ and 
ip' production in nucleus- nucleus collisions at the SPS ( A /s/v r /v ~ 17 GeV) and RHIC 
(\/SAriv=200 GeV) energies can be well described within the statistical hadronization 
model. We have recently shown [17] that there are crucial differences in elementary 
collisions compared to nucleus-nucleus for J/ip production. 

Besides the thermal parameters discussed above, which we keep unchanged, the 
model has as input parameter the charm production cross section in pp collisions, used 
to calculate the number of directly produced cc pairs N^ r which enter into the balance 
equation [H [15]: JV* = I^^jgSg} + g 2 c N% 

The centrality dependence of the nuclear modification factor R J Jj[ is shown in 
Fig. [3j The model reproduces very well the decreasing trend versus centrality seen 
in the RHIC data [18]. The larger R^t value at midrapidity is in our model due to the 
enhanced generation of charmonium around mid-rapidity, determined by the rapidity 
dependence of the charm production cross section. At the much higher LHC energy 
the larger charm production cross section could lead to a different trend as a function 
of centrality, depending on the magnitude of shadowing in Pb-Pb collisions. A generic 
prediction of the model is that the R 3 Jj± value at LHC is larger than at RHIC and 
this is confirmed by the preliminary ALICE data [19] measured at forward rapidity, 
which demonstrate, in our view, that charmonium is produced at LHC at the phase 
boundary (chemical freeze-out). If further confirmed by data (importantly, also on ip' 
production), this picture will consolidate the role of charmonium as a special observable 
to probe deconfinement of heavy quarks and to delineate the phase boundary of QCD 
matter with hadrons carrying heavy quarks. At the LHC, first results of T production 
have appeared [21\. The model predicts for T a suppression-like pattern of a similar 
magnitude as that of J/ip at RHIC. We predict p2] for Pb-Pb much smaller ratios 
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Figure 3. Centrality dependence of Rj/a f° r RHIC and LHC energies at mid-rapidity 
(left panel) and forward rapidity (right panel). The two curves shown for the LHC 
energy coorespond to a range of expected shadowing. The ALICE data shown in the 
right panel are preliminary results shown at this conference |19j . 



T(2S')/T(1S')=0.033 and T(25)/T(15')=0.005, compared to the values in pp collisions 
at Tevatron [20], 0.32 and 0.15, respectively. This feature is indicated in the CMS data 

PH. 
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